Screening model of metallic non-ideal contacts at integer quantized Hall regime 
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In this work, we calculate the electron and the current density distributions both at the edges 
and the bulk of a two dimensional electron system, focusing on ideal and non-ideal contacts. A 
three dimensional Poisson equation is solved self-consistently to obtain the potential profile in the 
absence of an external magnetic field considering a Hall bar defined both by gates (contacts) and 
etching (lateral confinement). In the presence of a perpendicular magnetic field, we obtain the spatial 
distribution of the incompressible strips, taking into account the electron-electron interactions within 
the Thomas-Fermi approximation. Using a local version of Ohms law, together with a relevant 
conductivity model, we also calculate the current distribution. We observe that the incompressible 
strips can reside either on the edge or at the bulk depending on the field strength. Our numerical 
results show that, due to a density poor region just in front of the contacts, the incompressible strips 
do not penetrate to the injection region when considering non-ideal contact configuration. Such a 
non-ideal contact is in strong contrast with the conventional edge channel pictures, hence has a strong 
influence on transport. We also take into account heating effects in a phenomenological manner and 
propose a current injection mechanism from the compressible regions to the incompressible regions. 
The model presented here perfectly agrees with the local probe experiments all together with the 
formation of hot-spots. 

PACS numbers: 73.23.Ad, 73.43.Fj, 73.43.-f 



I. INTRODUCTION 

Regardless of what system one is interested in theoret- 
ically 1 4 , the metallic contacts deposited on the two di- 
mensional charge systems are the most important ingre- 
dients of the measurement, both experimental!}^^ and 
theoreticall y 10 " 13 ^, where by contacts we mean both the 
injection/collector and probe contacts. In a standard 
Hall experiment (not necessarily quantized) an external 
magnetic field B is applied to the system in the direction 
normal of the plane. The excess electrons are injected 
from the source contact and are collected at the drain 
contact by applying a finite potential difference between 
these two contacts. Meanwhile, one measures the (elec- 
trochemical) potential difference either between the con- 
tacts along the sample (longitudinal potential difference 
Vu in the case of a fix imposed current /, resistance, Ri 
measured between A1-A2 or B1-B2 in Fig. [I]) or across 
the sample, the Hall potential Vjj (or resistance Rh, 
measured between Al-Bl or A2-B2 in Fig. [I]). Deter- 
mining the actual contact resistances of each element is 
a cumbersome matter. Experimentally, the usual way of 
obtaining the (non-ideal) contact resistance is the trans- 
mission line method, which gives an average approximate 
value by measuring a set of contacts defined on a stripe of 
mesa containing 2DES 14 . However, in the quantized Hall 
regime the determination of contact resistances becomes 
somewhat easier, since one a priori knows that the lon- 
gitudinal resistance should vanish and Hall resistance is 
quantized to an integer sub-multiple of the von Klitzing 
constant R K = h/e 2 = 25,812.807449(86) where e is 



the elementary charge and h is the Planck constant!^. 
The term non-ideal contact stands for a case where the 
edge-channels do not equilibrate at the contacts, hence 
reflection (or transmission) from channel to contact (or 
vice versa) is non-zero (not unity). In contrary, in an 
ideal contact all the channels are in equilibrium with the 
contacts, and reflection is zero. We should also clarify 
the notation ideal Ohmic contact: This term stands for 
a contact with finite resistance Rc, however, still obey 
the Ohm's law V = IRc where the applied current is 
/. Hence, a non-ideal contact for the Landauer-Biittiker 
formalism, i.e. non-zero reflection, can still be an ideal 
Ohmic contact. A very comprehensive review is provided 
in RefP. 

It is a formidable task to model the sudden change 
of the density of states (DOS) near the contacts, while 
at the metallic region the DOS is approximately infinite, 
whereas at the 2D system the DOS is either a constant 
(without magnetic field B, Dq) or v aries b etween no DOS 
at the Fermi energy (incompressible^^) and high DOS 
(compressible, with a degeneracy of eB/h), in the pres- 
ence of a high magnetic field. Recently, there were nu- 
merical efforts to model the contacts by graphical meth- 
ods^ considering interacting classical electrons or by the 
non-equilibrium network modeP. The ideal contacts are 
modeled as equipotential surfaces, where no density poor 
region resides just in front of the contacts. Whereas, the 
non-ideal contacts were described by different injection 
probabilities at the injection region by Kramer et a/, how- 
ever, density fluctuations near the transition region, as 
indicated by experiments^! which influences the current 
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distribution drastically is left unresolved. 

Here, we present our self-consistent ly calculated results 
which takes into account Landau quantization and di- 
rect Coulomb interactions at a mean-field approximation 
levepEU. 

In this paper, we investigate theoretically the current 
and charge density distributions considering quantizing 
high B fields and intermediate electron density concen- 
trations in the close proximity of the contacts, together 
with the entire sample. Our results point that, the as- 
sumption of ideal contacts cannot be justified, when con- 
sidering density inhomogeneities. As a guiding estimate, 
we employ the relevant formula provided by Chklovskii 
et al 20 to calculate the incompressible width just in front 
of the contact, 
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where, a\ is the width of first incompressible strip, a* B 
(~ 9.81 nm for GaAs/AlGaAs) is the effective Bohr ra- 
dius, Id is the depletion length, the bulk filling factor 
and \p is the Fermi wavelength. For a typical sample, 
i.e. Xf ~ 30 nm, and for v Q = 1.1 one observes that 
the depletion length should be less than 20-30 nm (4-5 
times the Bohr radius), which is much smaller compared 
to the experimental findings^. Hen ce, models based on 
ideal contacts has to be revisecPHS^ where the bulk elec- 
tron density cannot be reach within few a* B . Moreover, 
we present clearly the formation of hot-spots^ depend- 
ing on the current amplitude and field direction. We also 
take into account Joule heating following the lines of H. 
Akera and his co-workers, phenomenologically and pro- 
pose a self-consistent mechanism to inject current from 
compressible to incompressible regions (and vice versa). 
Our treatment of the contact regions is essentially based 
on the experimental findings of J. Weiss and his co- 
workers, namely scanning force microscopy experiments 
by E. Ahlswede^ and comprehensive structural investi- 
gation by O. Goktas^. Up to our knowledge, no com- 
prehensive study of contacts is accessible in the litera- 
ture that both take into account interactions, device ge- 
ometry, potential fluctuations and formation incompress- 
ible and compressible regions, which exist due to Landau 
quantization. 



II. THE MODELING OF CONTACTS 

In spite of the importance of the contacts at 2DESs, 
it was only recent that systematic investigations revealed 
the actual material composition at the close vicinity of 
the contacts^. They found that, the contacts annealed 
with the Ni/Ge/Au composition has the following prop- 
erties: i) The 2DES falls apart from the contacts at the 
Au-rich regions with an electron depleted region that ex- 
tends 50 to 150 nm. ii) In contrast, at the close proximity 
of NiGe-rich grains the 2DES can reach up to the metal- 
lic region, however with a density gradient till the bulk 




FIG. 1: (Color online) Schematic presentation of the Hall 
bar and contacts on a 3D GaAs/AlGaAs (red and blue re- 
gions) heterostructure (a). The sample has the dimensions of 
1.5 x 1.5 x 0.85 /im 3 , mapped on a matrix of 96 x 96 x 110 
embedded in a dielectric material, where open boundary con- 
ditions are assumed on the border of the dielectric material. 
Two donor layers, distributed homogeneously on the xy plane, 
provide electrons both to surface and to the 2DES. Contacts 
are modeled by i) A metallic gate on the surface ii) A metal- 
lic region inside the crystal iii) Fingers at the plane of 2DES, 
where fingers are to simulate bad contacts, (b) A sketch of 
2D projection of the Hall bar from top: Dark (gray and light 
gray) regions depict the contact, whereas the lateral confine- 
ment is shown by the etched regions. Color gradient repre- 
sents the electron density, together with Landauer-Biittiker 
edge states, in the case of an ideal contact (straight lines) 
and an non- ideal contact (broken curves). 



density n e \ is reached. The question how the current is in- 
jected to this region is answered by a model deduced from 
temperature dependent resistance experiments pointing 
that a very high and thin Schottky-barrier is formed at 
the contact/GaAs interface, hence the injection is due to 
tunneling. The most striking result of such an experimen- 
tal finding is that the Biittiker type edge-states cannot 
equilibrate within the contacts considering an electron 
poor region just in front of the contacts. For sure, for thin 
(Schottky) barriers and no poor region LB formalism is 
adequate and can be extended to thick barrier regime. 
However, the electron density gradient in front leads to 
almost full reflection of the Biittiker type edge-states, as 
depicted in Fig. [TJd. 

In the above work it was pointed that a model which 
also takes into account such a density gradient in front 
of the contacts is admired, especially when considering 
high perpendicular magnetic fields, where compressible 
and incompressible strips form. Here, we present such 
a model by solving the Poisson equation in 3D self- 
consistently to obtain potential and charge density distri- 
butions near the contact, together with the entire sample, 
employing a forth order grid techniqu e 18 * 24 ! In the next 
step we use the potential profile obtained at the plane of 
2DES as an initial condition to calculate same quantities 
in the presence of a perpendicular B field. Our calcula- 
tion scheme is based on a mean-field Thomas Fermi ap- 
proximation improved by spatial course-graining to sim- 
ulate quantum mechanical effects, such as the finite ex- 
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tent of the wavefunctiondU At a last step we impose 
a fixed current in y direction and calculate the current 
distribution utilizing a local version of the Ohm's law^. 
Here, we assume that the density of states (DOS) and the 
temperature of electrons do not vary considerably on the 
quantum mechanical length scales, i.e. these quantities 
are position independent. The conductivity model is ob- 
tained from the self-consistent Born approximation^!! 
(SCBA) and local conductances are assumed to be di- 
rectly related with the local electron density. Such an 
approximation is somewhat crude, however, is valid if the 
charge density changes slowly on the correlation length 
of the remote impurities 26 . An improved calculation 
scheme^, already supports our assumption on the con- 
ductivity model. Moreover, our results are independent 
of the choice of the conductivity model, since we are not 
interested in the details of plateau to plateau transition 
regions. The only necessary ingredient of our conduc- 
tivity model is to have vanishing (in fact exponentially 
small, at zero temperature) longitudinal conductivity at 
incompressible regions and a Hall conductance propor- 
tional to the electron density. The last requirement also 
implies that, the Hall conductance is quantized at the 
incompressible regions. 

In Fig. we show schematic presentation of the sys- 
tem at hand. The 2DES lies some 150 nm below the sur- 
face, where surface potential is pinned to the mid gap of 
GaAs. Narrow Hall bar is defined by etching at the sides 
by a depth similar to 65 nm, whereas the metallic fingers 
reside at the plane of 2DES. To simulate the effects result- 
ing from annealing, we also introduced a metallic region 
65 nm below the surface. The contact region is high- 
lighted in Fig. [I}}, where dark regions correspond to gold 
alloy and light regions represent NiGe alloy. The front 
of gold alloy is strongly depleted, shown by the white 
region at the 2DES, whereas, an electron poor stripe re- 
sides next to the NiGe alloy. For further references, we 
define the finger widths (wf) and lengths (If) differently 
to simulate good (wf — nm and If = nm) and bad 
contacts (wf > 65 nm and If > 190 nm). In the case 
of good contacts we assume that the depleted regions in 
front of the Au grains is negligibly small, whereas for the 
bad contacts the depleted region is generated with the 
metallic fingers by keeping the potential at -0.45 V and 
there are no fingers at the regions corresponding to NiGe 
grains. The strong confinement in x direction is modified 
by the contact region. We begin our discussion with a 
case which has an one to one correspondence in classical 
Hall effect: We set the B field such that only the lowest 
Landau level is partially occupied, hence, behaves like 
a metal. Throughout this paper, following the previous 
works, we neglect the spin degeneracy, since the essential 
physics is independent of the origin of the single particle 
gap, however, spin gener alized versions of the screening 
theory already existP^SI. It is common to describe the 
2D electronic system in the presence of a perpendicular 
field by a dimensionless parameter; the filling factor, 
which measures the occupancy of the Landau levels below 
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FIG. 2: (Color online) Spatial distributions of the filling fac- 
tor v(x,y) (a), together with the current density (b) as a 
function of lateral coordinates. Color scale denotes density 
gradient, whereas arrows (blue) present the amplitude and 
direction of the imposed excess current with an amplitude 
< 0.2 nA, guaranteeing liner response regime. The calcula- 
tions are performed at Tiuj c /Ef = 4.38 and default tempera- 
ture kT/EF ~ 0.08. The electron poor region is a result of 
metallic contacts at the surface and penetrated metal, kept 
at -0.5 V. The associated electrochemical potential distribu- 
tions are shown in the inset of (b) , following the vertical solid 
and broken lines, respectively. Boxes indicate the hot-spot 
regions, separated from the contact regions by broken (green) 
vertical lines. Here, contact region is defined by the density 
gradient, hence Landau levels also do form in this region. The 
current is injected exactly from y = /xm line. 



the Fermi energy given by the ratio of electron number 
density n e i to the magnetic flux density n^. One can ex- 
press the filling factor also in terms of magnetic length 
Ib (= yJfi/eB) as v — 27rl B n e \ and the local version as 
v(x,y) = 2iti B riQ\(x, y). Therefore in this case the filling 
factor is below 2 (due to spin degeneracy) and there are 
many available states at the Fermi energy, similar to a 
metal. Such a case (also locally) is called as compressible, 
where screening is nearly perfect and electrons can be re- 
distributed according to applied electric fields. Fig. [2^i 
presents the filling factor distribution (or electron den- 
sity) calculated for the geometry described above, how- 
ever, without metallic fingers at the plane of 2DES to 
mimic an ideal Ohmic contact. Due to the strong con- 
finement by the etching, a large electron depleted region 
resides on both sides of the Hall bar. Note once more 
that, such a contact with a density poor region is not an 
ideal contact considering Landauer-Biittiker edge chan- 
nel picture, as depicted in Fig. (TJd. One observes that, 
the electronic density is slightly small in front the con- 
tacts compared to the bulk. If we consider the classi- 
cal Hall effect, the current (and electrochemical poten- 
tial) distribution is well knowrPHl there will be two spots 
where electrochemical potential is highest (depending on 
the field direction) at the right-bottom and lowest at the 
left-top corners. For an illustrative demonstration we 
suggest the reader to check Fig.l of Ref.[T2land also Fig. 4 
to compare with actual first principle many-body simu- 
lations, performed at the classical regime together with 
interactions. The experimental findings also support the 
formation of these spots^^^. Here the current ampli- 
tude is measured in units of evpkp j (j\L y \ where vp, 
are the Fermi velocity and wave vector, respectively, and 
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L y (= 1500 nm) is the length of the sample, hence the 
current amplitude is less than 0.2 nA directed in positive 
y direction. The current distribution also presents two 
regions where most of the transport takes place, denoted 
by (red) boxes in Fig. [2}}, we observe that the current 
is almost homogeneously injected by the source contact, 
however, is confined to these regions. Once the current 
enters to the bulk of the sample, it is distributed homoge- 
neously following th e loca l electron density, likewise the 
classical Drude resultP^M] This regime overlaps perfectly 
with the local probe experiments 6 , where we also show 
the demonstrations of the electrochemical potential dis- 
tribution across the sample in the lower inset of Fig. (2}d. 

Next, we lower the B field such that a wide incom- 
pressible region (denoted by black) is formed at the bulk 
of the sample, Fig. |3^i. We also show the corresponding 
current distribution in Fig. [3Jd, where a positive current 
is imposed in y direction. The hot-spots are still visible, 
however, now the current is injected to the incompressible 
bulk instead of a compressible bulk. Injecting current to 
an incompressible region is somewhat counter intuitive 
at a first superficial look. Before elucidating this, we 
would like to emphasize that our results presented are 
calculated self-consistently for a given background po- 
tential and imposed current, i.e. finite electrochemical 
potential difference between the source and drain con- 
tacts, together with the local Ohm's law. Moreover, we 
do not assume that the current should be injected to the 
incompressible regions. However, self-consistency results 
in such a case. In fact the key concept is the electro- 
chemical potential difference between contacts: We inject 
electrons which have an electrochemical potential energy 
larger than the Fermi energy otherwise (namely, if they 
were all at the Fermi energy) there would be no net cur- 
rent. It is useful to re-introduce 3 the definition of (local) 
Ohm's law; 



V/x*(r)/e = p(r)j(r)=E(r), 



(2) 



where /i* (r) is the position dependent electrochemical po- 
tential, p(r) is a two by two tensor describing the local 
resistivities and j(r) is the local current density together 
with the local electric field E(r). The calculation proce- 
dure is such: First we obtain the electrostatic properties 
of the system at total equilibrium, (i.e. /i* is position in- 
dependent and equals to Ep), then impose a small (com- 
pared to equilibrium) electrochemical potential difference 
between source and drain contacts and re-calculate the 
electron distribution depending on newly calculated to- 
tal potential. Once local equilibrium conditions are sat- 
isfied together with the numerical convergence, namely 
the potential and electron density distributions do not 
change within the numerical accuracy of 10 _T in the last 
iteration, the imposed electrochemical potential differ- 
ence is increased step by step to the target value. In 
fact, in our calculations presented here, we imposed a 
very small current amplitude therefore, the density dis- 
tribution does not change substantially. In other words, 
we are in the linear response regime, however, we are 
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FIG. 3: (Color online) Spatial distribution of the electron 
densities (left panel) and current densities, same as Fig. [2] 
calculated at selected B field strengths, Tiuj c /Ef — 3.78 (a- 
b), 3.02 (c-d) and 2.26 (e-f) at default temperature, (e-f) 
mimics an ideal contact in edge-channel picture, where no hot- 
spots are observed. Insets at the right panel demonstrate the 
electrochemical potential distributions similar to the previous 
figure. 



not limited to this regime and our result s in the non- 
linear regime are published elsewher e 1 18 * 33 *. Moreover, if 
one takes into account heating effect^ due to the high 
current density at the hot-spots, it is easy to see that 
the incompressible strips at the injection and emission 
regions are melted, hence the excess current can be trans- 
ferred between compressible and incompressible regions 
(and vice versa) easily. Including such local tempera- 
ture effects are far complicated to be incorporated within 
our calculation scheme. However, one can still estimate 
the local electron temperature variation at the hot-spots 
from the work by H. Akera, namely Ref. [2j focusing on 
Fig.l and Fig. 2b (as a function of filling factor and lat- 
tice temperature) or as a function of imposed current 
amplitude (Fig. 7). Depending on the parameters we 
employed here, we estimate that, the electron tempera- 
ture at the hot-spots are 5-10 (and even more in strip 
case) percent higher than the lattice temperature, which 
essentially implies that the incompressible region at the 
hot-spot melts, i.e the single particle energy gap closes 
due to large derivative of the Fermi function. In the in- 
set of Fig. (3}d, we also depict the spatial variation of the 
associated electrochemical potential. It is observed that, 
at the hot spot region ji*(x) varies in a highly nonlin- 
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FIG. 4: (Color online) Self-consistently calculated carrier den- 
sity (a) and current density distributions (b), where (b) fo- 
cuses on the 2DES region. The fluctuations at the contour 
plot present the anomalies near the contacts, smearing out 
the quantized Hall effect. 



ear manner, having a maximal variation exactly at the 
spot (solid thick line), whereas, far from the injection 
contact potential presents an s-shape behavior (broken 
line) similar to the experimental 5 and theoretical 32 find- 
ings. Hence, the Hall resistance measured at the center 
is quantized, however, near the hot-spots deviate from 
the plateau value considerably. Fig. [3^ and Fig. [3jl, cor- 
responds to a magnetic field strength where a deformed 
incompressible ring is formed close to the edges. Noth- 
ing different occurs at injector (and collector) region at 
this B value, the non-linear potential drop is observed 
with slight modifications. In contrast the electrochemi- 
cal potential profile changes substantially at the center, 
the strong variation only takes place where the incom- 
pressible strips reside, inset Fig. [3^. We observe that, 
once the current is injected to the incompressible strip it 
is kept being confined to this region due to the absence 
of back-scattering, till the next hot-spot. One deduces 
from the electrochemical potential profile that the Rh 
is quantized, similar to our previous works^ and experi- 
ments 5 . 

So far we have observed the features of an ideal Ohmic 
contact (non-ideal in the language of Landauer-Biittiker) 
due to the electron poor region in the front of contacts. 
However, our results show that the current can be di- 
rectly injected to the incompressible strips through the 
compressible region at the contacts. We have seen that, 
near the contact regions the essential futures of the quan- 
tized Hall effect is lost. In contrast far from the contacts 
it is recovered, we think that these numerical results are 
in accord with the ones of Ref. |13| and with the experi- 
mental findings 5 9 . Now the question is, what happens if 
an incompressible edge strip touches the contact region, 
i.e an ideal contact. Such a case occurs as a natural re- 
sult of the self-consistency, which we present in Fig. [3^. 
The direct Coulomb interaction and Landau quantiza- 
tion (and Lorentz force), pins the electron density just in 
front of the contact to an (even) integer filling factor at 
2.26huj c / Ep. The corresponding, current distribution is 
shown in Fig.[3j, where one cannot observe the formation 
of the hot-spots. Since, all the current is directed along 
the incompressible strip starting from the injection con- 
tact (namely, due to the absence of back-scattering). In 



this case, one can think of an ideal contact in the termi- 
nology of Landauer-Biittiker edge-state formalism. Un- 
fortunately, such a case is in contrast to the local probe 
experiments^ and other experiments which report hot- 
spotsP. 



The influence of the contact quality on the quantized 
Hall effect was also investigated in Ref. [9j It was re- 
ported that, measurements performed on deeper lying 
heterostructures considering contacts which are defined 
in bad contacting direction, the plateaus are not well de- 
veloped. In this last part we model bad contacts, yet 
utilizing thick metal fingers residing in the plane of the 
2DES, namely we set the finger length to be 190 nm 
and width 65 nm, at the injection and collector channels. 
These metallic fingers essentially represent electronically 
depleted regions in front of Au grains. Whereas, the 
Ni/Ge grains correspond to regions in between the fin- 
gers (kept at ground), where a density gradient occurs. A 
typical result is shown in Fig.[4| considering same param- 
eters chosen in Fig. [3^ and Fig. [3j, together with contour 
lines of the current density. Here, we superimposed den- 
sity fluctuations induced by negatively charged (-0.5 V) 
metallic fingers to simulate effects resulting from large Au 
grains, both in the longitudinal (corresponding to injec- 
tion and collector contacts) and lateral (simulating probe 
contacts) directions, where latter is defined by thicker 
fingers (If = 530 nm and Wf = 80 nm). The first obser- 
vation is that, the current is no longer injected (almost) 
homogeneously from the contact to the 2DES, moreover, 
no incompressible strips form at the inner parts of the 
fingers due to the strong potential variation. Hence, scat- 
tering is enhanced. To simulate the effect of bad side con- 
tacts we also impose metallic fingers on the sides. We set 
the widths and lengths similar to the injection/collector 
contacts, whereas in Fig. [2] and Fig. [3] these lengths were 
set to zero. We observe that, the density in the lat- 
eral direction also presents oscillations and the widths 
of the incompressible strips vary considerably, even they 
become narrower than the extent of the wavefunction. 
Hence, the back-scattering free strips are lost, resulting 
in non-quantized Hall plateaus. Instead of hot-spots, we 
see that the current is diverted to the right-bottom cor- 
ner due to Lorentz force. The influence of bad probe con- 
tacts would be suppressed if these side contacts reside far 
apart. Unfortunately, our numerical abilities are limited 
by the dimension of the matrix mapping the Hall bar, 
therefore we cannot perform systematic investigations in 
this regime. In any case, our results suggest that, once 
the density and potential fluctuations imposed by probe 
contacts are suppressed at large distances the ideal be- 
havior would be recovered. Another important point to 
note is, once the sample becomes wider and the effect of 
potential fluctuations are suppressed by screening near 
the side contacts, the quantized Hall effect would be re- 
covered, since one only measures the (electrochemical) 
potential difference between the side contacts. 
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III. CONCLUSIONS 

Here, we investigated the influences of contacts on the 
current and density distributions, considering a Hall bar 
under quantized Hall conditions, within the screening 
theory. We modeled the contacts following the ideas put 
forward by O. Goktas based on the experimental find- 
ings, namely the formation of an electron poor region 
in the close proximity of contacts. It is shown that, in 
the classical regime where all the system is compress- 
ible and behaves like a metal, the well known hot-spots 
are formed at the corners of the sample. At lower fields, 
where Landau quantization becomes important, incom- 
pressible strips (or regions) are formed and current is 
injected to these states via hot-spots in the case of ideal 
Ohmic contacts. Here we also discussed heating effects in 
a phenomenological manner and argued that due to high 
current densities at the hot-spots the incompressible re- 
gions melt locally. Moreover, we showed that due to the 
electron poor region, Landauer-Biittiker edge-channels 
are reflected from the contacts and can be considered 
as non-ideal contacts. We also provide numerical results 
such that if the B field and the contact structure is chosen 
appropriately, one can still obtain an ideal contact in the 
terminology of Landauer-Biittiker formalism. However, 
one cannot observe the formation of hot-spots in this 
case which is in contradiction with the experiments men- 
tioned, unfortunately. In a final discussion, we also simu- 
lated bad contacts taking into account potential fluctua- 
tions resulting from Au/Ni/Ge alloys by placing metallic 
fingers in the plane of the 2DES. We observed that, the 
incompressible strips are destroyed and large amount of 
scattering takes place both near the injection/collector 



and side probe contacts, hence the quantized Hall ef- 
fect is lost. Our findings, we think, are in accord with 
the experimental findings and also with the numerical 
investigations considering the classical regime and non- 
equilibrium network model. The investigation of actual 
injection of electrons through the Schottky barrier and 
local temperature effects demands a more complicated 
calculation than presented in this work, however, our nu- 
merical investigations support the idea that such an in- 
vestigation might modify the picture presented slightly 
in a small parameter space. 

While submitting our manuscript, we have encoun- 
tered a very recent experimental paper reporting on the 
anisotropic depletion at contact interfaces measured by 
scanning force microscope^. The reported potential pro- 
files strongly support our model, which identifies pres- 
ence of an incompressible strip in front of ohmic con- 
tacts, that might decouple compressible bulk from the 
contacts. Moreover, our modeling of good/bad contacts 
also perfectly agree with their findings. 
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